*This is the replication do file for:

******* Knutsen, Carl Henrik (2010). Investigating the Lee thesis: How bad is democracy for Asian economies? ****************************** 
************************************ European Political Science Review. 2(3): 451-473. ****************************************************  


clear
**Insert own location for data file in the below expression:
*use "M:\Aførsteam\replication\EPSR_repl.dta", clear

keep country AggregFHI GDPPCgrowthWDI grownrollmenttertiary grownrollmentsecondary grownrollmentprimary grossinvestmentpercgdp WDIGDPPP logGDPPPP asia politysquared _stack year latam  dec85 dec87 dec91 dec92 Politynonanarchy africass	AREG	asiapac	budhistkonftaoshinto	budhkonfutaoshinto	catholic dec82	dec83	dec84	dec86	dec88	dec89	dec90	dec93	dec94	dec95	dec96	dec97	dec98	dec99long	easteursov	hindu	indigeneous islamshia	islamsunni	lead2growthgdpppc	lead2infergrowth logdur	logdurpl1	loggdppp	loggdpppinclinfer	loginferpop	logpop	logWDIpop	mideanafr	orthodox	orthodoxarm	pluspolity	protestanglic	protestanglicmethod	romcath	shia	sunni	time18701913	time19141945	time19461972	time19722003	westeurplusoffshots	wgicorrupt

**This operation could be performed to shave the dataset down to only Asian countries. However, note that some analyses below are conducted on global samples (with interaction terms for Asian countries)
*gen includ = 0
*replace includ = 1 if asiapac==1
*replace includ = 1 if asia==1
*keep if includ==1
*drop includ

label variable lead2growthgdpppc "GDP per capita growth"
label variable logpop "Log of population size"
label variable logdur "Log of the political regime’s duration"
label variable logWDIpop "Log of population size (WDI)"
label variable grossinvestmentpercgdp "Gross investment per gdp"
label variable grownrollmenttertiary "Enrolment tertiary"
label variable grownrollmentsecondary "Enrolment secondary"
label variable grownrollmentprimary "Enrolment primary"
label variable logGDPPPP "log GDP per capita (WDI)"
label variable GDPPCgrowthWDI "GDP growth per capita (WDI)"
label variable loggdppp "log of per capita GDP (Maddison)"
label variable Politynonanarchy "Polity coefficient"
label variable sunni "Sunni Muslim dummy"
label variable hindu "Hindu dummy"
label variable shia "Siha dummy"
label variable catholic "Catholic dummy"
label variable protestanglic "Protestant"
label variable orthodox "Orthodox"
label variable budhkonfutaoshinto "Buddhism, Confucian, Taoist, Shinto"
label variable indigeneous "Indigenous religion"
label variable time18701913 "time dummy 1870-1913"
label variable time19141945 "time dummy 1914-1945"
label variable time19461972 "time dummy 1946-1972"
label variable time19722003 "time dummy 1973-2003"
label variable AggregFHI "Freedom House Index"
label variable AREG "ACLP democracy measure"
label variable loggdpppinclinfer "log of per capita GDP (Maddison) including interpolated data"
label variable loginferpop "Log of population size including interpolated data"
label variable logdurpl1 "Log of (political regime's duration +1)"
label variable country "Country name"

**NEWDATASET
generate leadinvestpergdp = grossinvestmentpercgdp[_n+2]
generate leadprimschool = grownrollmentprimary[_n+2]
generate leadsecschool = grownrollmentsecondary[_n+2]
generate leadtertschool = grownrollmenttertiary[_n+2]
generate leadGDPPCgrowthWDI = GDPPCgrowthWDI[_n+2]

set more off

* table 1 p. 464
xtpcse lead2growthgdpppc logpop logdur sunni hindu budhkonfutaoshinto Politynonanarchy loggdppp time18701913 time19141945 time19461972 time19722003 if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
estimates store OLSPCSEperiod
xtpcse lead2growthgdpppc logpop logdur sunni hindu budhkonfutaoshinto Politynonanarchy loggdppp dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
estimates store OLSPCSEdecade
xtpcse lead2infergrowth loginferpop logdur sunni hindu budhkonfutaoshinto Politynonanarchy loggdpppinclinfer time18701913 time19141945 time19461972 time19722003 if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
estimates store OLSPCSEperiodlarge
xtpcse lead2infergrowth loginferpop logdur sunni hindu budhkonfutaoshinto Politynonanarchy loggdpppinclinfer dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
estimates store OLSPCSEdecadelarge
estout OLSPCSEdecade OLSPCSEperiod OLSPCSEdecadelarge  OLSPCSEperiodlarge, cells(b(star fmt(%9.3f)) t(par fmt(%9.2f))) starlevels (^ 0.10 * 0.05 ** 0.01 *** 0.001) stats (N)style (tex)


* table 2 p. 465
**Post-45 sample**
xtpcse lead2growthgdpppc logpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdppp dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1 & year>1945, correlation(psar1) rhotype(tscorr) pairwise
estimates store a

**Post 60
xtpcse lead2growthgdpppc logpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdppp dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1 & year>1960, correlation(psar1) rhotype(tscorr) pairwise
estimates store b
estout a b, cells(b(star fmt(%9.3f)) t(par fmt(%9.2f))) starlevels (^ 0.10 * 0.05 ** 0.01 *** 0.001) stats (N)style (tex)

* table 3 p. 465
**FE
xtreg lead2growthgdpppc logpop logdur Politynonanarchy loggdppp time18701913 time19141945 time19461972 time19722003 if asiapac==1, fe robust
estimates store FEperiod
xtreg lead2growthgdpppc logpop logdur Politynonanarchy loggdppp dec82 dec83 dec84 dec85 dec87 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1, fe robust
estimates store FEdec
xtreg lead2infergrowth loginferpop logdur Politynonanarchy loggdpppinclinfer time18701913 time19141945 time19461972 time19722003 if asiapac==1, fe robust
estimates store FEperiodlarge
xtreg lead2infergrowth loginferpop logdur Politynonanarchy loggdpppinclinfer dec82 dec83 dec84 dec85 dec87 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1, fe robust
estimates store FEdeclarge
estout FEdec FEperiod FEdeclarge FEperiodlarge, cells(b(star fmt(%9.3f)) t(par fmt(%9.2f))) starlevels (^ 0.10 * 0.05 ** 0.01 *** 0.001) stats (N)style (tex)

estout OLSPCSEdecade OLSPCSEdecadelarge FEdec FEdeclarge, cells(b(star fmt(%9.3f)) t(par fmt(%9.2f))) starlevels (^ 0.10 * 0.05 ** 0.01 *** 0.001) stats (N)style (fixed)

* table 4 p. 466
**Post45&60
xtreg lead2growthgdpppc logpop logdur Politynonanarchy loggdppp dec82 dec83 dec84 dec85 dec87 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1 & year>1945, fe robust
estimates store a
xtreg lead2growthgdpppc logpop logdur Politynonanarchy loggdppp dec82 dec83 dec84 dec85 dec87 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1 & year>1960, fe robust
estimates store b
estout a b, cells(b(star fmt(%9.3f)) t(par fmt(%9.2f))) starlevels (^ 0.10 * 0.05 ** 0.01 *** 0.001) stats (N)style (tex)

generate Dec50s = 1 if year>1949& year<1960
generate Dec60s = 1 if year>1959& year<1970
generate Dec70s = 1 if year>1969& year<1980
generate Dec80s = 1 if year>1979& year<1990
generate Dec90s = 1 if year>1989& year<2000
generate Dec00s = 1 if year>1999
mvencode Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s, mv(0) override

*WDI p.466 pcse
xtpcse leadGDPPCgrowthWDI islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig Politynonanarchy logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise
*FHI
xtpcse leadGDPPCgrowthWDI islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AggregFHI logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise
*AREG
xtpcse leadGDPPCgrowthWDI islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AREG logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise

*WDI p.466 FE
xtreg leadGDPPCgrowthWDI islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig Politynonanarchy logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, fe robust
*FHI
xtreg leadGDPPCgrowthWDI islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AggregFHI logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, fe robust
*ARGE
xtreg leadGDPPCgrowthWDI islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AREG logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, fe robust

generate asiapolityinteract = asiapac*Politynonanarchy // generating interaction Polity and Asian regional dummy

**interaction global sample
xtpcse lead2growthgdpppc logpop logdur westeurplusoffshots easteursov africass asiapac mideanafr latam sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdppp time18701913 time19141945 time19461972 time19722003 asiapolityinteract, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2growthgdpppc logpop logdur westeurplusoffshots easteursov africass asiapac mideanafr latam sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdppp dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long asiapolityinteract, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2infergrowth loginferpop logdur westeurplusoffshots easteursov africass asiapac mideanafr latam sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdpppinclinfer time18701913 time19141945 time19461972 time19722003 asiapolityinteract, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2infergrowth loginferpop logdur westeurplusoffshots easteursov africass asiapac mideanafr latam sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdpppinclinfer dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long asiapolityinteract, correlation(psar1) rhotype(tscorr) pairwise

* p. 467 interaction term between polity and log GDP 
generate interactionpolloggdp =loggdppp*Politynonanarchy
label variable interactionpolloggdp "interaction between polity and log GDP"
generate interactionpolloggdpinfer =loggdpppinclinfer*Politynonanarchy
label variable interactionpolloggdpinfer "interaction between polity and log GDP" //sjekk
xtpcse lead2growthgdpppc logpop logdur interactionpolloggdp westeurplusoffshots easteursov africass asiapac mideanafr latam sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdppp dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2infergrowth loginferpop interactionpolloggdpinfer  logdur westeurplusoffshots easteursov africass asiapac mideanafr latam sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdpppinclinfer dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long, correlation(psar1) rhotype(tscorr) pairwise

**Chow test p. 467 ec growth 
generate post1980 = dec98+dec99long
generate chowpost80polity = post1980*Politynonanarchy
generate chowpost80FHI = post1980*AggregFHI
generate chowpost80AREG = post1980*AREG

generate chowpost90polity = dec99long*Politynonanarchy
generate chowpost90FHI = dec99long*AggregFHI
generate chowpost90AREG = dec99long*AREG

xtpcse lead2growthgdpppc logpop chowpost80polity logdur sunni hindu budhkonfutaoshinto Politynonanarchy loggdppp time18701913 time19141945 time19461972 time19722003 if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2growthgdpppc logpop chowpost80polity logdur sunni hindu budhkonfutaoshinto Politynonanarchy loggdppp dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2infergrowth loginferpop chowpost80polity  logdur sunni hindu budhkonfutaoshinto Politynonanarchy loggdpppinclinfer time18701913 time19141945 time19461972 time19722003 if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2infergrowth loginferpop chowpost80polity logdur sunni hindu budhkonfutaoshinto Politynonanarchy loggdpppinclinfer dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise

**Investment analyses p. 468
xtpcse leadinvestpergdp islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig Politynonanarchy logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse leadinvestpergdp islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AggregFHI logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse leadinvestpergdp islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AREG logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise

xtreg leadinvestpergdp islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig Politynonanarchy logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, fe robust
xtreg leadinvestpergdp islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AggregFHI logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, fe robust
xtreg leadinvestpergdp islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AREG logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, fe robust

*Chow-test
xtpcse leadinvestpergdp chowpost80polity islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig Politynonanarchy logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse leadinvestpergdp chowpost80FHI islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AggregFHI logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse leadinvestpergdp chowpost80AREG islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AREG logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise

**Testing Barro hyp
xtpcse lead2growthgdpppc politysquared logpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous pluspolity loggdppp time18701913 time19141945 time19461972 time19722003 if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2growthgdpppc politysquared logpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous pluspolity loggdppp dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2infergrowth politysquared loginferpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous pluspolity loggdpppinclinfer time18701913 time19141945 time19461972 time19722003 if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2infergrowth politysquared loginferpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous pluspolity loggdpppinclinfer dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise

xtpcse lead2growthgdpppc politysquared logpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous pluspolity loggdppp dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2infergrowth politysquared loginferpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous pluspolity loggdpppinclinfer dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long, correlation(psar1) rhotype(tscorr) pairwise

**Testing interaction
xtpcse lead2growthgdpppc interactionpolloggdp logpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdppp time18701913 time19141945 time19461972 time19722003 if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2growthgdpppc interactionpolloggdp logpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdppp dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2infergrowth interactionpolloggdpinfer loginferpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdpppinclinfer time18701913 time19141945 time19461972 time19722003 if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse lead2infergrowth interactionpolloggdpinfer loginferpop logdur sunni shia catholic protestanglic orthodox hindu budhkonfutaoshinto indigeneous Politynonanarchy loggdpppinclinfer dec82 dec83 dec84 dec85 dec86 dec87 dec88 dec89 dec90 dec91 dec92 dec93 dec94 dec95 dec96 dec97 dec98 dec99long if asiapac==1, correlation(psar1) rhotype(tscorr) pairwise

*Human capital p.469
xtpcse leadprimschool islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig Politynonanarchy logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse leadprimschool islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AggregFHI logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise
xtpcse leadprimschool islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AREG logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr) pairwise

xtpcse leadsecschool islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig Politynonanarchy logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr)hetonly pairwise
xtpcse leadsecschool islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AggregFHI logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr)hetonly pairwise
xtpcse leadsecschool islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AREG logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr)hetonly pairwise

xtpcse leadtertschool islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig Politynonanarchy logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr)hetonly pairwise
xtpcse leadtertschool islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AggregFHI logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr)hetonly pairwise
xtpcse leadtertschool islamsunni islamshia romcath protestanglicmethod orthodoxarm hindu budhistkonftaoshinto indig AREG logGDPPPP logdurpl1 logWDIpop Dec50s Dec60s Dec70s Dec80s Dec90s Dec00s if asia==1, correlation(psar1) rhotype(tscorr)hetonly pairwise

*Insert own location for data file
**save "M:\Aførsteam\replication\EPSR_repl.dta", replace   